Diamagnetism versus Paramagnetism in charged spin-1 Bose gases 
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It has been suggested that either diamagnetism or paramagnetism of Bose gases, due to the charge 
or spin degrees of freedom respectively, appears solely to be extraordinarily strong. We investigate 
magnetic properties of charged spin-1 Bose gases in external magnetic field, focusing on the competi- 
tion between the diamagnetism and paramagnetism, using the Lande-factor g of particles to evaluate 
the strength of paramagnetic effect. We propose that a gas with g < 1/V8 exhibits diamagnetism 
at all temperatures, while a gas with g > 1/2 always exhibits paramagnetism. Moreover, a gas with 
the Lande-factor in between shows a shift from paramagnetism to diamagnetism as the temperature 
decreases. The paramagnetic and diamagnetic contributions to the total magnetization density are 
also calculated in order to demonstrate some details of the competition. 
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I. INTRODUCTION 

The Bose gas plays a significant role in understanding 
a series of exotic quantum phenomena, including super- 
fluidity and superconductivity. It is well known that the 
ideal Bose gas exhibits the Bose-Einstein condensation 
(BEC) below a critical temperature. In 1938, London 
first connected BEC to the A transition in Helium.— In 
1946, Ogg proposed that the BEC of bosonic electron- 
pairs might result in superconductivity. 2 Furthermore, 
Schafroth 3 - and Blatt and Butler- showed that an ideal 
gas of charged bosons exhibits the essential equilibrium 
features of a superconductor. Although the BCS theory^ 
revealed that electrons in a superconductor form Cooper- 
pairs, as opposed to real-space pairs, the Schafroth-Blatt- 
Butler theory is helpful to the understanding of supercon- 
ductivity and more importantly it stimulated research 
interest in charged Bose gases. 

The charged Bose gas (CBG) is solely of academic in- 
terest in its own right. Especially, it exhibits nontriv- 
ial magnetic properties. For example, the condensation 
phenomenon in CBGs is strongly affected by the exter- 
nal magnetic field, owing to the quantization of the or- 
bital motion of charged particles in magnetic neld£>£~— 
Schafroth clearly indicated that an arbitrarily small value 
of the magnetic field introduces qualitative changes: BEC 
does no longer occur, 3 May extended this idea to a d- 
dimensional CBG and found that it can condense only for 
d > 4. 7 This point was reexamined by other researchers 
based on different methods! 10 ! 11 Meanwhile, the orbital 
motion results in extremely large Landau diamagnetism 
in CBGs. It was pointed out that the 3-dimensional CBG 
displays Meissner effect at low temperatures i 3 - ! 12 ' 13 More 
recently, Alexandrov quantitatively accounted for the en- 
hanced normal-state diamagnetism of superconducting 
cuprates using the CBG models 

On the other hand, neutral bosons with spin (spinor 
bosons) have also been studied theoretically^ The re- 
alization of spinor Bose condensate in optically trapped 
alkali atoms 16 stimulates new research interest. The con- 
stituent atoms, such as 87 Rb, 23 Na, and 7 Li have (hy- 



perfine) spin degrees of freedom and thus a magnetic 
moment. Their spin degrees of freedom become active 
in purely optical traps and thus investigation of their 
magnetic properties becomes possible. Yamadai^, and 
Simkin and Cohenii calculated the magnetization of neu- 
tral spinor bosons in magnetic field and found that once 
BEC takes place, the magnetization remains finite even 
if H = 0, as if the system was magnetized spontaneously. 
The zero-field susceptibility tends to diverge as the tem- 
perature goes down to the BEC temperature. Moreover, 
rigorous proofs presented by Eisenberg and Lieb show 
that the magnetization and zero-field susceptibility at fi- 
nite temperatures exceed that of a pure paramagnet. 18 
All the results indicate that the neutral spinor bosons 
take on extraordinary paramagnetic effect in magnetic 
field. 



Bearing in mind the two contrary features of the Bose 
gas, one immediate question which must be raised is: 
which kind of magnetism will manifest itself if the bosons 
possess both the charge and spin degrees of freedom? 
Analogously, this question has been answered for charged 
Fermi gases, e.g., the electron gas, where both the para- 
magnetism and diamagnetism are relatively weak. For 
the electron gas, the diamagnetic part of the zero-field 
susceptibility is one-third (in absolute value) of the para- 
magnetic part, so altogether the gas is paramagnetic. For 
charged spinor Bose gas, since both the paramagnetism 
and diamagnetism can be extremely large, their compe- 
tition is a more interesting problem. 



In this paper, we study the competition between para- 
magnetism and diamagnetism of a charged spin-1 Bose 
gas in external magnetic field. In Section II, a model con- 
sisting of both the Landau diamagnetic effect and Pauli 
paramagnetic effect is proposed. The total magnetiza- 
tion density as well as its paramagnetic and diamagnetic 
parts are calculated respectively. Section III presents a 
detailed discussion of the obtained results. In Section IV, 
a brief summary is given. 
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II. THE MODEL 

The orbital motion of a charged boson with charge q 
and mass m* in a constant magnetic field B is quantized 
into the Landau levels, 
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where j = 0,1,2,... labels different Laudau levels and 
lu = qB/(m*c) is the gyromagnetic frequency. Since 
uj cc B, u) can be used to indicate the magnitude of the 
magnetic field in the following discussions. We assume 
the magnetic field is in the z direction. Each Landau 
level is degenerate with degeneracy equal to 



A 



qBL x L y 
2Tthc 



(2) 



Here we suppose the gas is in a box with Li — ¥ oo, where 
Li is the length of the box in the ith direction. For a spin- 
1 boson, the Zeeman energy levels split in the magnetic 
field due to the intrinsic magnetic moment associated 
with the spin degree of freedom, 
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where a refers to the spin-z index of Zeeman state 
\F = l,m_p = a) (<7 = +1,0,-1) and g is the Lande- 
factor. The quantization of the orbital motion and the 
Zeeman effect give rise to the Landau diamagnetism and 
Pauli paramagnetism, respectively. 

We consider an assembly of N bosons, whose effective 
Hamiltonian reads 
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where fi is the chemical potential. The charged spinor 
bosons have been discussed theoretically in the context 
of relativistic pair creation^. However, magnetism of 
charged spinor Bose gases is less studied and of major 
interest in the present work. 

The grand thermodynamic potential is formally ex- 
pressed as 
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where /3 = (fcsT) . Converting the sum over k z to 
continuum integral, we get 
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where V is the volume of the system. Then using the 
Taylor expansion and performing the integral over k z , 
we have 
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This treatment has been used by Standen and Toms^i 
to deal with scalar Bose gases. As they mentioned, this 
theory is more reliable at high temperature. Similarly, we 
introduce some compact notation for the class of sums, 



la /2 e -lx(e+S) 

£ ^ tt ^HL (l-e-'*)" ' 



(8) 



where x = f3hui and fi + ghqoB j (to* c) = (^ — e)huj. With 
this notation we may rewrite Eq. ([7]) as 
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where D — 3. Then the density of particles n = N/V 
can be derived from the thermodynamic potential, 
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The magnetization density M is written as 
1 fdQ 
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It is convenient to introduce some dimcnsionless param- 
eters, such as M = m*cM/(nhq), uj = fuv/iksT*), 
t = T/T* and x' — ul/t, to re-express equations <jTU]) 
and (JTTJ, 
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FIG. 1: (a) The total magnetization density (M), (b) the 
paramagnetization density (M p ), and (c) the diamagnetiza- 
tion density (Md) as a function of g for fixed magnetic fields 
at t — 0.1. Here the solid line, dashed line and dotted line 
correspond to Zo — 0.05, 0.3, and 3, respectively. Inset: m as 
a function of g. 



Here the characteristic temperature T* is given by 
k B T* = 2irh 2 n 2 / 3 /m* . The Bose-Einstein condensa- 
tion temperature of spin-1 Bose gas with density n is 
just defined as k B T c = 27rft 2 n 2 / 3 /{m*[3C(3/2)] 2 / 3 } « 
/csT*/3.945. The dimensionless notation EJJar, S] should 
be 

00 la/2 -lx'(e+8) 

km=y: (1 i e - lX ' )K . ( i4 ) 

z=i 1 ' 

where // = fi/(k B T*) and /z' + gcruJ = (| — e)oJ. The 
two variables // and M are determined by Eqs. (fT2|) and 

(USD. 

III. RESULTS AND DISCUSSIONS 

In our model, both the charge and spin degrees of free- 
dom are taken into account, which are described respec- 
tively by the Landau energy term e^ fc and the Zeeman 
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FIG. 2: Plots of the critical value of Lande-factor, g c as a 
function of 1/t for fixed values of Zo. The field is chosen as: ZJ 
= 10 (dash-dot-dotted line), 3 (dotted line), 0.5 (short dotted 
line), 0.3 (dashed line), 0.1 (short-dashed line), and 0.05 (solid 
line). 

energy term e z a e in the Hamiltonian (Q| . In the case that 
the Lande-factor g tends to zero, the model degenerates 
into a charged scalar boson model which exhibits strong 
diamagnetism as already intensively discussed^— As g 
becomes larger, the paramagnetic effect is strengthened. 

We calculate the dimensionless magnetization density 
M as a function of g, as shown in Fig. Ufa). M is neg- 
ative in the small g region, which means that the dia- 
magnetism dominates. The absolute value of M is larger 
in the stronger field uJ. For each given value of ZJ, M 
increases monotonically with g. M changes its sign from 
negative to positive at a critical value of g, noted as g c 
hereinafter, reflecting that the paramagnetism becomes 
dominant as g increases. Note that the slope of the M 
curve is dependent on g. When g is near to zero, M in- 
creases slowly but the slope rises quickly with g. It means 
that the interplay between diamagnetism and paramag- 
netism is complex and nonlinear. However, in the strong 
paramagnetic region, M grows almost linearly with g. 

Figure |Tfb) plots the paramagnetic contribution to M 
(named as the paramagnetization density), M p = gm, 
with an inset showing m — n\ — and Figure [TJc) 
shows the diamagnetic contribution to M (named as the 
diamagnetization density), = M — M p . The increas- 
ing tendency of M p is similar to that of M. It is notewor- 
thy that the diamagnetization density is not suppressed, 
but enhanced as g becomes larger. Comparing Figs[TJa), 
[Tfb) andQJc), it can be seen that the increase in M comes 
from the paramagnetic effect. In the small g region, both 
M p and Md arc strengthened nonlincarly with increas- 
ing g. As shown in the inset of Fig. [ljb), m grows very 
quickly. Nevertheless, both m and Md flatten out in the 
large g region. So the slope of M curve is mainly due to 
the paramagnetization. 

According to discussions above, the critical value of the 
Lande-factor, g c , is an important parameter to describe 
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FIG. 3: Curves of g c — 1/t obtained respectively from the 
Bose-Einstein (BE) and Maxwell-Boltzmann (MB) statistics 
with fixed values of £7. Where U = 10 (dash-dot-dotted line, 
BE; solid line, MB), 3 ( dotted line, BE; short dotted line, 
MB), and 0.3 (dashed line, BE; short dashed line, MB). 



the competition between the diamagnetism and param- 
agnetism. g c is a function of the temperature t and the 
magnetic field 57. Figure [5] shows g c as a function of 1/t 
in different magnetic fields 57 = 10, 3, 0.5, 0.3, 0.1 and 
0.05. Obviously, g c varies monotonically with the tem- 
perature t, while its dependence on the field 57 is not 
simple. For example, in the low temperature region, g c 
decreases with decreasing 5J at a given temperature as 57 
is still larger than 0.3, then it rises up as 57 goes down 
further from 57 0.3. 

As already mentioned, the results of our theory are 
more credible at high temperature. In the high temper- 
ature limit, g c seems universal with respect to different 
choices of magnetic field, g c \t-^oo ~ 0.35356. For a given 
magnetic field, g c increases as the temperature falls down. 
This suggests that the diamagnetic region is larger at low 
temperatures than at high temperatures. Although the 
exact value of g c can not be obtained at very low temper- 
atures, its variation trend can be estimated from Fig. [2] 
It seems that g c ranges from 0.475 to 0.50 in various 
magnetic fields. 

It is useful to reexamine the high temperature behav- 
iors of g c by generalizing above calculations to a spin-1 
Boltzmann gas, since the Bose-Einstein statistics reduces 
to Maxwell-Boltzmann statistics in the high temperature 
limit. The grand thermodynamic potential based on the 
Maxwell-Boltzmann statistics reads 
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Then equations of the dimensionless chemical potential // 
and magnetization density M are derived respectively, 
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FIG. 4: Shown are plots of dimensionless M as a function of 
1/t for each given value of uJ at a fixed g. The value of ZJ 
were in sequence as 10 (dash-dot-dotted line), 3 (dotted line), 
0.3 (dashed line), and 0.05 (solid line), (a) corresponds to 
g — 0.35. (b) corresponds to g = 0.45. (c) corresponds to 
ff = 0.5. 



and 
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Substituting Eq. JTS]) into Eq. ([I7|. yields 
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An analytical formula for g c can be obtained, 
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The value of g c can be derived from Eq. (TK)1) exactly 
in two limit cases: <? c |t->oo = l/v8 ~ 0.35355 and 
<?c|t-K) = 1/2- The value of g c for a Boltzmann gas is 
reasonably equal to that of a Bose gas in the high tem- 
perature limit. Figure |3] plots the numerical solutions of 
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FIG. 5: The dimensionless M p and M<j as a function of 1/t for fixed values of w. From left to right g = 0.35, 0.45 and 0.5, 
respectively. For each given value of g, the field is chosen as: ZJ = 10 (dash-dot-dotted line), 3 (dotted line), 0.3 (dashed line), 
and 0.05 (solid line), respectively. 



Eq. (fTT))) and compares with the Bose gas. An interesting 
point is that the low-temperature-limit value of g c is also 
consistent with that of a Bose gas at low temperature but 
in high magnetic field, although the Maxwell-Boltzmann 
statistics is just valid in the high temperature region. 
The stronger the field, the better the accordance. 

The g c — 1/t curves in Figs. [2] and |3] mark the bound- 
ary between the diamagnetic and paramagnetic regions. 
The gas exhibits diamagnetism at all temperatures and 
in all magnetic field when g < g c |t^oo ~ 0.35355, while 
always paramagnetism when g > 0.5. Whereas, the mag- 
netic properties seem more complicated in the interme- 
diate region, g c \t^oo < g < 0.5. Figures SJa-c) illustrate 
the dimensionless magnetization density M for the three 
different cases, respectively. 

Figure a) denotes the case of g = 0.35 < g c \t^oo- 
The temperature-dependence of M is very similar to that 
of charged scalar Bose gases,— as if the paramagnetic ef- 
fect associated with the spin degree of freedom is thor- 
oughly hidden. The diamagnetism is even stronger at 
lower temperatures. As the external field tends to be 
weak, a sharp bend appears gradually on the curve, which 
located at the point corresponding to the BEC tempera- 
ture in zero field^ In our model, the BEC temperature 
for a spin-1 gas is (l/t) c ~ 3.945. Figure |4|c) shows M 
for a paramagnetic case to the contrary when g = 0.5. 
M is always positive in the field at all temperatures. An 
interesting phenomenon is that the M — 1/t curve shows 
up a peak in this case. The decline in M at low tem- 
peratures is attributed to the diamagnetic effect. When 
weakening the magnetic field, the peak is lowered and 
moves to low temperatures. Fig. Hlb) depicts the case 



with g in the intermediate region, which looks quite sim- 
ilar to Fig.[4}c). The key difference is that M can change 
its sign from positive to negative as the temperature de- 
creases, indicating that the system undergoes a shift from 
paramagnetism to diamagnetism. 

Figures HI demonstrate the total magnetic performance 
of the charged spin-1 Bose gas. We now turn to exam- 
ine the underlying paramagnetic and diamagnetic effects 
for each corresponding case. As shown in Figs. [5j both 
the paramagnetization density and the diamagnetization 
density become more stronger at lower temperatures. As 
the external field is reduced, the strengthening of M p 
and Md becomes fast near the temperature close to the 
BEC point in zero field. As g grows, Md is only slightly 
strengthened but M p increases significantly and thus it 
can go beyond Md. M p thoroughly exceeds Md when 
g rises to 0.5. If g increases further, the paramagnetic 
effect can be so strong as to cover up the diamagnetic 
effect completely. Then the total magnetization density 
M becomes monotonously increasing with lowering the 
temperature and finally reaches a plateau, instead of a 
peak, at low temperatures. 

IV. SUMMARY 

This paper has studied the interplay between param- 
agnetism and diamagnetism of the ideal charged spin-1 
Bose gas. The Lande-factor g is introduced to describe 
the strength of paramagnetic effect caused by the spin 
degree of freedom. The gas exhibits a shift from diamag- 
netism to paramagnetism as g increases. The critical 
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value of g, g c , is determined by evaluating the dimen- 
sionless magnetization density M, Our results show that 
g c increases monotonically as t decreases. In the high 
temperature limit, g c goes to a universal value in all dif- 
ferent magnetic fields, g c \t^oo = 1/V8- At low temper- 
atures, our results indicate that g c ranges from 0.475 to 
0.50 as the magnetic field varies. Therefore, a gas with 
g < l/v8 exhibits diamagnetism at all temperatures, 
but a gas with g > 1/2 always exhibits paramagnetism. 
In cases where l/y/8 < g < 1/2, the Bose gas under- 
goes a shift from paramagnetism to diamagnetism as the 
temperature decreases. 

In order to depict some details of the competition 
between paramagnetism and diamagnetism, the para- 
magnetic and diamagnetic contributions to the total 
magnetization density are also calculated. No doubt 
that the paramagnetism is enhanced with increasing g. 
Surprisingly, the diamagnetism is not suppressed, but 
slightly strengthen. This implies that the competition 
between para- and dia-magnetism is nontrivial. When g 
is fixed, both the paramagnetism and diamagnetism be- 
come stronger as t decreases. As in the scalar case, there 
is no Bose-Einstein condensation when the magnetic field 
is present, no matter how small it is. However, evidence 
of the condensation can be seen in the magnetization 
density as the magnetic field is reduced. 

At last, we briefly discuss experimental aspects possi- 



bly relevant to the present work. Although the charged 
spinor Bose gas has not been realized so far, the up- 
to-date achievement in experiments makes it attainable 
perhaps in the near future^ For example, it is already 
possible to create ultracold plasmas by photoionization of 
laser-cooled neutron atomsi^i The ions can be regarded 
as charged bosons if their spin is an integer. The Lande- 
factor for different magnetic ions could be different. As 
reported by Killian et al. 21 , the temperatures of elec- 
trons and ions are as low as 100 mK and 10 /zK, respec- 
tively. One can expect that the temperature could be 
lowered near to the BEC temperature with new advance- 
ment in experimental techniques. The diamagnetism- 
paramagnetism effects manifest themselves near the BEC 
point. On the other hand, several ferromagnetic super- 
conductors have been discovered since 2000*22, Cooper 
pairs in such materials are likely in the spin-triplet state, 
thus behave somewhat like charged spin-1 bosons. Note 
that the charged spin-1 boson model can not be applied 
directly to describe triplet superconductors, but it could 
offer some help in understanding magnetic properties of 
such materials. 
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